Spin decoherence due to Fluctuating Fields 



Piotr Szarikowski, 1,2 M. Trippcnbach, 1 and Y. B. Band 2 

'institute of Theoretical Physics, University of Warsaw, ul. Hoza 69, PL-00-681 Warszawa, Poland 
2 Department of Chemistry, Department of Physics and Department of Electro- Optics, 
and the Use Katz Center for Nano-Science, Ben-Gurion University, Beer-Sheva 84105, Israel 

The dynamics of a spin in the presence of a deterministic and a fluctuating magnetic field is 
solved for analytically to obtain the averaged value of the spin as a function of time for various 
kinds of fluctuations (noise). Specifically, analytic results are obtained for the time dependence of 
the expectation value of the spin, averaged over fluctuations, for Gaussian white noise and Guassian 
colored noise, as well as non-Gaussian telegraph noise. Fluctuations cause the decay of the average 
spin vector (decoherence). For noise with finite temporal correlation time, a deterministic component 
of the field can suppress decoherence of the spin component along the field. Hence, decoherence can 
be manipulated by controlling the deterministic magnetic field. A simple universal physical picture 
emerges which explains the mechanism of the suppression of the decay. 
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I. INTRODUCTION 

A quantum system is almost invariably coupled to an environment that degrades its coherence. Decoherence of 
quantum systems is a fundamental problem, with implications across all branches of physics. The temporal evolution 
of a spin system, in the presence of coupling to an environment is perhaps the simplest phenomenon involving 
decoherence. For example, a particle of spin 1/2 [which can be mapped onto any two- level system (TLS) [l|, and 
vice versa] has served as a useful paradigm [2 J6j. Specific examples include the decoherence of a qubit in a quantum- 
computer or quantum-information processor J decoherence of an electron spin in ESR or a nuclear spin in NMR 0] , 
spin dynamics in a magnetometer [9(, accuracy and stability limitations of atomic clocks due to interactions with an 
environment fioj , decoherence of two- level quantum dot systems [TP . [HJ , and nitrogen vacancy centers in diamond 
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The main goal of this manuscript is to discuss the decoherence of a spin system caused by an environment. Under 
certain assumptions, which we discuss in the following section, the interaction of a spin with an environment can 
be represented by evolving the spin in an effective magnetic field, £?( off ) = £? + Bsit). Here Bq is a deterministic 
magnetic field (which can be time-dependent) that can be used to exert control over the spin. The magnetic field 
BE{t) models the influence of an environment on the system, and is represented by a vector stochastic process b(t). 
Averaging over fluctuations corresponds to tracing out the environmental degrees of freedom. This yields a reduced, 
nonunitary dynamics wherein the averaged spin decoheres in time. 

The physical properties of the environment determine the statistical properties of -Be(£), which in turn determine 
the type of stochastic process b(t). Here we consider several types of stochastic fluctuations. A prototype model 
for fluctuations is Gaussian white noise [HI, [l9| , wherein the random process has vanishing correlation time. As a 
consequence, the system driven by it evolves without memory, since it has no way to probe the past. A more general 
type of noise is Gaussian colored noise [2(J Hl[, wherein the random process has a finite correlation time, hence, the 
system will have memory. A particular case of colored noise that is not Gaussian, telegraph noise, has been studied 
in connection with quantum dot qubit systems A main result of our analysis is that, when the spin evolves 

with memory, it can be manipulated by applying a sufficiently strong deterministic field Bq to significantly suppress 
decoherence. 

The outline of the this paper is as follows. Section [TTl develops the equations of motion for a spin in the presence 
of a deterministic magnetic field and a stochastically fluctuating magnetic field. Section UTT1 describes evolution of the 
system with Gaussian colored noise fluctuations. In that section we present a simple, intuitive view of the nature 
of the decay of the averaged spin due to field fluctuations. We also show that this decay can be suppressed when 
the deterministic magnetic field is strong enough, and the correlation time is long enough (i.e., when the system has 
memory). Section ITVl develops the white noise limit of colored noise. Section fVTl considers telegraph noise, which is a 
specific from of non-Gaussian noise. Finally, a summary and a conclusion are presented in Sec. I VIII 
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II. SPIN DYNAMICS IN A FLUCTUATING MAGNETIC FIELD 



We consider fluctuations in an otherwise deterministic system that results from the application of a effective random 
field generated by the environment. The back-coupling of the spin to the source of the noise is not taken into account 
(Van Kampen [18( calls this external noise). Hence, the stochastic properties of the noise result only from the 
environment which is unaffected by the system. We assume that these properties can be measured or otherwise 
deduced. For each particular system, one should carefully check whether the back-action can indeed be neglected. 
Examples of systems for which back-action was neglected include: nitrogen- vacancy centers in diamond affected by 
magnetic moments of impurity spins [l5l Il6j , electron spin in a quantum dot affected by nuclear magnetic moments 
that fluctuate due to crystal lattice vibrations [l2l Il3j . and magnetic noise in atom chips caused by fluctuations of 
electron currents in wires that make up the atom chip [22| . 

Consider a particle of spin S with magnetic moment /2 = fiS in the presence of a time-dependent field B(t). The 
dynamics is determined using a Zeeman-typc Hamiltonian 

H = —ft ■ B{t) = &(t) ■ S, (1) 

where S is the vector of spin operators satisfying commutation relations [Si,Sj] = iheijkSk, f2(£) = —/j,B(t)/h is the 
Rabi frequency vector, B(t) is a magnetic field, fx = g/j,o is the magnetic moment, g is the g-factor and fiQ is the Bohr 
(or nuclear) magneton. The magnetic field B{t) is the sum of a deterministic field Bo, whose direction defines the 
z-axis, and a fluctuating field b(t), hence the Rabi frequency is given by 



n(t) = + Uy {t) , (2) 




where Slo = —fiBo/h and u)(t) = —fib(t)/h. Wc have taken the deterministic field to have only a component along the 
z-axis for simplicity. 

In the Heisenberg representation, the evolution of the spin operators is given by the equations, 



^ H \t) = -[H(t),Sr>(t)\, (3) 
where the superscript (H) denotes the Heisenberg picture. In our case, the Hamiltonian contains a stochastic term, 

H(t)=n Q (t)S z +uj(t)-S. (4) 

The commutator on the right hand side of the Heisenberg equation of motion, Eq. is easy to evaluate, and leads 
to the vector form of the equation, 

^-S {H \t) = {Q (t)e z +uj(t)} x S^(t), (5) 



1 



AH) 



dt 



and the initial condition is simply 




(6) 



For convenience, let us take the deterministic field Qq to be constant in time. Then, the solution to Eq. ([5]) with 
the initial condition ((6]) can be written as a 3x3 evolution matrix, ?7(t,0), acting on the initial vector, S^ H \0). It is 

convenient to write the equations, 5j H) (t) = Y%=i U lJ {t,Q)S i J H) (0), in vector notation as S^(t) = U(t,0)S^ H \0), 
where the hat " indicates a 3x3 matrix. Substituting this form into the equations of motion yields 

d U(t 7 0)= (j((,0). (7) 



dt 



fe=i 



Here the matrix elements of i h are defined as = ejjy, where e^jy is the Levi-Civita symbol. The matrices e* satisfy 
commutation relations similar to those of angular momentum operators, 

[i i ,i j ) = e ijk £ k . (8) 



3 



Consequently, they are generators of rotations in a real, three-dimensional vector space. For example, e z generates 
rotations around z axis, 



R(8,z) ee exp{6>e z } = sin(9 cos(9 . (9) 




The evolution operator U is a rotation operator that causes precession of the spin vector around the instantaneous 
rotation axis defined by the total field, and the instantanous Rabi frequency vector is f2(t) = f^o +t3(t). Hence U is a 
stochastic matrix. The equation of motion of the spin does not depend upon the details of the representation of the 
spin, hence the solution to the dynamics for an arbitrary spin is also given by (J7J . The only indication of the spin 
representation, i.e., the dimension 25 + 1 of the representation, is in the initial condition. 

Experimental measurements of the spin will of necessity correspond to quantum expectation values and averages over 

the stochastic fluctuations of the magnetic field. Hence the quantity (S(t)) corresponds to experimental measurements 
of the spin. Here, the symbol { ) means quantum average, and the symbol means average over stochastic variables. 
Note that the average over fluctuations and the average over the initial quantum state factorizes in the following way: 

(W)) = U(t,0)(S(0)). 

For spin 1/2 (i.e., for a TLS), the average expectation value of the spin S = ha/2 fully determines the state of 
the system. This is because the density matrix g can be expressed using the Pauli matrices a which, together with 
identity matrix 1, form a basis in the space of 2x2 hermitian matrices. In the Schrodinger representation, 

^=i(l + Tr{^^(0)}-<?) =\(l+Wfl-Z)- (10) 

and the expectation values of the spin at time t is given by vector parameter A(i) = (<j(t)). Thus, at any time t, there 
are three parameters, Ai(i), i = x,y,z, which completely specify the density matrix. This does not mean that the 
dynamics of the density matrix for any (2S + l)-level system can be represented in terms of the dynamics of a particle 
with spin S coupled to a magnetic field, since the number of parameters necessary to specify the density matrix larger 
than the three components of the spin [24| . 



III. GAUSSIAN COLORED NOISE 

Initially, we consider a process that is Gaussian, isotropic and Markovian The field fluctuates in all three 

dimensions with independent components (isotropy) and statistical properties of each component is completely de- 
termined by the first two moments (Gaussian), the average and the correlation function. The stationary condition 
means that the correlation function depends on the time difference and the average is time independent. According 
to Doob's theorem [26|, if the process is also Markovian, the correlation function is exponential and the average is 
zero, i.e., 

^(i) = 0, (11) 

Wi(i)wj(t') = n(t - t') 5ij ^ule-^ 6 i:i . (12) 

Here the bar indicates average over the fluctuations, n(t — t') is the correlation function with correlation time r c , the 
average fluctuating Rabi frequency vanishes, uj = 0, and the variance of the Rabi frequency, Wg, specifies the strength 
of the fluctuations. 

Several comments regarding Guassian colored noise are in order. First, note that the Gaussian process can model 
a wide variety of environments due to the central limit theorem; if the environment consists of many individual and 
independent fluctuating elements, the central limit theorem substantiates the approximation of using a Gaussian 
process as the effective field affecting the spin. Secondly, if there is no back-action of the spin on the environment, 
and the environment is in a stationary state (e.g., is in equilibrium) then the process is stationary as well. Finally, the 
isotropy assumption regarding the random magnetic field produced by the environment is appropriate if no particular 
special direction can be associated with the environment. 

To solve the equation of motion ([7]), first we eliminate the so called drift term due to the deterministic part of the 
magnetic field, fio e z U{t, 0), from the right hand side of the equation by transforming to the rotating frame. Then 
Eq. ([7]) takes the form 



!#'(*,()) = (*(*)• !>(*)) u'(t,o). 



(13) 
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In this equation we introduced transformed matrices, denoted by a prime, e'(t) and U'(t, 0), which are given by 

(e%(t) = Y, R im{t^,z)e k mn R nj {-m ,z) = R k i{m ,z)e\ j , (14) 

Ui J {t ll t 2 ) = ^ RimjhQo, z) U mn (ti, h) Rnji-hftp, z), (15) 

where R(ttto, z) = exp{tilo£ z } is a deterministic Euclidean rotation matrix that results in counterclockwise rotation 
around the z-axis by the angle fi * [see Eq. ©]. Note that the scalar product oj(t) ■ i'(t) in Eq. ([13)) can be rewritten 

as 

u{t) ■ £'{t) = (R(-t Q , z) iu(t)j ■ 1 = Lj'{t) ■ i. (16) 

The time-dependent rotation matrix R(ttto,z) "mixes" field fluctuations in the x and y component. In the Lab 
frame, u> x and lo v arc independent, however, because of the mixing, they are correlated in the rotating frame and the 
correlation functions are modified in the following way: 



4(*iH(* 2 ) = -«i(tiK(t 2 ) = «(ti - h) sin [n Q (h - 1 2 )\ , (17) 



u/Mw'M = W'(tiK(*2) = K(tl-t2)cOB[no(tl-ta)]. (is) 



wi(*iM(*a) = «(ti - t a ) <S:k. (19) 

The stochastic field in the rotating frame, u:'(t), remains a Gaussian process, since it is a linear combination of Gaussian 
processes. Note that any time- independent rotation does not change the correlation functions of the isotropic Gaussian 
vector process uj(t). 

The formal solution to Eq. (|13[) is given in the terms of the time-ordered exponential function, 

U'(t,0) =rexp|^e 1 ^ <4(r)drj , (20) 

defined by the power series, Texp J Q A(t)cIt = 1+ L cItiA(ti) + L dr\ /J" 1 dr 2 A(ti)A(t 2 ) + . . ., where T is the time- 
ordering operator (32j . The average value of L/', t/'(t, 0), can be computed by noting that it has the form of a moment 
generating functional {Mx\ £(t)] \ for matrix processes, Xj(i) = u[{t)i l , i = x,y,z, where £(f) = (6r(i), £,y(t), ^ z (t)) is 
the vector of trial functions [38|. Then, the average can be "absorbed" inside the exponential function by expressing 
Mx [£(£)] through a cumulant generating functional {Kx[£,{t)\\ in terms of cumulant averages denoted here by a 
double bar (see Ref. [25| for a discussion of the cumulant expansion methods used here) : 



^{r)dr\ = MxUt)\ 



U'(t,0)=Texp^2 it I 

= r exp{K ue [i]} = r cx P jrcxp f o ^ 

= Tcxp J Ve l e 3 f d 

{ i,3 J ° 



(r)dr 



*-2W^(ri)w;(T 2 )j. (21) 

In the last step, we took advantage of the fact that all cumulants beyond the second vanish for Gaussian processes, 
and u^(t)a^(i') = u;£(t)u;j(t') if the average vanishes. 

Substituting the correlation functions given by Eqs. (|17|) . (|T51) and (fT5|) . we find the following expession for the 
cumulant generating functional K evaluated at £(i) = 1, 

K^l) = (e z ) 2 [ dr 7 (r) + ((£ x ) 2 + (e^f) [ dr 7+ (r) + {e x e» - f dr 7 _(r), (22) 
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where the rates 7 are given by integrals of the correlation functions: 



7+W 



dTUj' z (T)u;'M = / dTK(T), 

Jo 

t rt 



cLtlo' x {t)u:' x (0) 



dTU> y {T)uji y {0) 



dr k(t) cos(51o r ); 



7-(<) = / dTUj' x (T)oj'JO) = - dr W ;(rK(0)= / dr k(t) sm(n Q r) 



(23) 
(24) 
(25) 



K can be further simplificcl if we take into account properties of the rotation generators e % . From the commutation 
relations ((HJ), the commutator in the third term of Eq. (f22j) is equal to e z . It is easy to verify that the sum of squares 
of e x and e y in the second term of if can be written as i x2 +i y2 = — (21 + e 2 ), hence, the parts ofEq. (f22j) commute 
with each other. Thus, the average evolution matrix, which is equal to the exponential of K, factorizes into three 
parts, 



U'(t,0) =Tcxp{AV £ -[l]} = cxp<|V 2 / dr 7 (r) - (21 + e z2 ) f dr 7+ (r)+£ z f dr 

I Jo Jo Jo 



7_(r) 



e ~ r « 
e" r « 
1 



decay caused by 
the fluctuations in z 



e- r +W 
e- r +C) 
e - 2r +(*) 



cosT_(i) -sinr_(i) 
sinr_(i) cosT_(i) 
1 



(26) 



decay caused by 
the fluctuations in xy plane 



precession frequency modification 
due to fluctuations in xy plane 



where 

r+(t) = 

r_(t) = 

r(t) = 



dr7+(r) 



cJt7_(t) = 



dr 7 (r) = w 2 t c 

















(r 





i + n 



t — r n 



O i + 
l-e-*/' 
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2 T 2 



e T c cos riot 
t_ 

s T <= sin f2ot - 



2ft T, 



1 



2£1 T C 



"O'c 



t c sin Qq^ 
e~~ cos flot — 1 



(27) 
,(28) 
(29) 



The first term on the right hand side of Eq. ([26)) describes the decay caused by the fluctuations in the direction 
parallel to the applied field f2o = Qo&z- These fluctuations affect only spin components perpendicular to the z axis, 
and the decay rate T(t) does not depend on the intensity of the constant field. The second term originates from the 
fluctuations in the xy plane, which is perpendicular to the applied field. Fluctuations in the plane couple to all the 
components of the spin and cause their decay. The key feature is that the rate of decay due to these fluctuations is 
dependent on f2g- Finally, the third term, also related to the fluctuations in the xy plane, describes the modification 
of the precession frequency of the system caused by the noise. 
The averaged expectation values of the spin are given by: 



(S x (t)) = e-[ r W+ r +W]{cos[r!oi + r_(i)] sin [fiot + r_(t)]<S„)} > 

Wfi) = e-PW+r+^lsiniaot + r.^^^+cosIOot + r.^)]^)}, 
iMtj) = e~ 2r +V{S z ). 
In the limit O -» 0, we obtain T + (t) -4 T(t) and r_(i) -> 0, hence Eqs. ([30" |) -([32 )l reduce to 

- 2r W(S). 



(S(t)) 



f2n->0 



-> e 



(30) 
(31) 
(32) 

(33) 



Substituting Eq. (|29|) into (|33[) yields the analytic form ([33| of the isotropic decay of the average spin, regardless of 
the initial state of the spin, as can be seen in Fig. [1] For times t ;g> r c , the decay is simply exponential with the rate 
of decay given by T(t) « WgT c i. For times t r c , the rate of decay is F(t) « Wgt 2 /2. In the absence of an external 
field, the fluctuations in all three directions remain independent and the average evolution matrix decomposes in the 
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f/T, 



FIG. 1: The expectation value of the spin, (S(t)), subject to Gaussian colored noise, versus time for Qo — and ujqt c = 1. The 
analytic formula for (S(t)) is given in Eq. ()33[) . 




t/T c 



FIG. 2: (S z (t)) subject to Gaussian colored noise versus time for nonvanishing constant ilo (see legend) and ojqt c = 1. Equation 
(|32|l gives the analytic formula for (S z (t)). 



following way: 



£1 =0 




e- r W 



U(t,0) = e" r W 1 e~ r W . (34) 



o o e - r (') 




fluctuations in z fluctuations in y fluctuations in x 



Here the first part comes from fluctuations in z and causes the decay of the x and y components of the spin, the 
second part comes from fluctuations in y and causes the decay of the perpendicular components, and the third part 
originates from fluctuations in x direction causes the decay of the y and z components of the spin. 

Now let us consider finite H . The dependence of (S z (t)), the component parallel to the external field f2 , versus 
time for three values of dimensional parameter f2nr c is plotted in Fig. [2] As we discussed previously, z component of 
the average spin decays due to fluctuations in the xy plane, which are affected by the presence of the constant field. 
For ttoT c <C 1, the quantity (S z (t)) is clearly very similar to the results plotted in Fig.[U since r + (t) — > T(t). As £Iqt c 
increases, the decay of (S z (t)) is significantly slowed, and r + (t) — > in the limit of ftoT c —> oo, i.e., the decay is fully 
suppressed. 

Figure [3] plots (S x (t)) versus time for the case when (S y (0)) — and Fig. [4] presents the time dependence of the 

length of the vector (S±(t)) = ((S x (t)), (S y (t)),0). The decay of the components of the average spin perpendicular to 

the constant field fin, is not only induced by the fluctuations in the xy plane, but also by those in the z direction [see 
Eq. (Prj|) ]. We have previously seen that the effects of the fluctuations in x and y are inhibited if Oqt c is large enough. 
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L0 1 2 3 4 5 

the 

FIG. 3: The average value (S x ), when the spin is subject 
to Gaussian colored noise for nonvanishing constant Slo 
(see legend) and loot c = 1. Equation (|30l) together with 
initial condition (S y (0)) = yield the analytic formula for 
the function plotted. The average value (S y ) for initial 
condition (S x (0)) — is identical to the results shown 
here. 




the 



FIG. 4: Average value of j(oj_(t))[ = 

y {Sx(t)) + (S y (t)) subject to Gaussian colored 
noise for nonvanishing constant fio (see legend) and 
ojoTc = 1. Equations (|30[) and (|31l) yield the analytic 
formula for the function plotted. 



However, the influence of the fluctuations in the direction parallel to applied field is unperturbed by the presence of 
Oo. As long as the z component of the field is fluctuating the perpendicular component of the average spin will decay 
to zero for any value of £Iqt c as t/r c — > oo. For small OoT c the decay is, again, similar to the isotropic case shown 
in Fig. ([2]), but when Oot c increases, the decay is modulated by oscillations which result from the precession of the 
spin about the external field Oo with frequency modified by the fluctuations. For t ^> t c , the oscillation frequency is 
given by Oo ^1 + 1 j^pfp J [see Eq. (|28|)]. When Oot c 3> 1 the frequency of precession tends to Oo as the frequency 
modification due to fluctuations, r_(t) vanishes, in a similar manner as that of r + (i). 

For spin 1/2 system (a TLS), the decay of (S x (t)) and (S y (t)) correspond to the decay of the coherence of the 
system, whereas the decay of (S z (t)) corresponds to the vanishing of population imbalance. If the constant field is 
taken along x (as opposed to along z), decay of (S x (t)) is suppressed, and the decoherence time T2 ~ r c Oo/o;o, which 
can be made arbitrarily large (suppressed decay) by increasing Oq. 



IV. THE WHITE NOISE LIMIT 



The white noise limit of colored noise involves going to the regime of extremely rapid fluctuations of the field, i.e., 
vanishing correlation time r c . Unfortunately, simply taking r c — > does not yield the correct limit. If we put r c = 
in Eq. (|29p. we obtain that the decay rate of the average spin equals zero, and fluctuations do not affect the evolution 
of the spin. The proper limit can be obtained if the strength of fluctuations, as given by the variance of the field, 
Lo(t) 2 = uiq, is increased as the correlation time t c goes to zero, so that u>qT c remains constant. Thus, we arrive at the 
limit when the correlation function of the colored noise becomes proportional to a Dirac delta function, 

Wi(t)wj(*') = oj 2 e-^Sij = 2lo 2 t c (J- e'^ j % ^ T ^ D } 2 DS(t - t')5 lv (35) 

\Zt c J t c -s-0 

In the limit of k(t) — > 2DS(t) the decay rates T + (<) = T(t) = 2Dt and r_(t) = 0. The decay of the expectation value 
of the spins are then isotropic and arc not affected by the presence of a constant field Oo, i.e., 

/ e -r(*)-r+(t) C osr_(i) -e- r W- r +Wsinr_(t) \ 
(S(t))= e- r (*)- r +(*)sinr_(t) e- r W- r +(*)cosr_(t) i?(O t, z) (S) -> e- 4Dt R(n t, z) (S). (36) 

V e- 2r +(*) / 

Thus, the spin decay is isotropic, purely exponential with decay rate 4D, and is independent of Oo. Moreover, the 
precession due to Oo is unaffected by the fluctuations. In the white noise limit, the precession due to deterministic 
part of the field and the decoherence due to field fluctuations are completely decoupled. 
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V. MEMORY EFFECTS 



In the expression for (S(t)) for the case of Gaussian colored noise, the decay rates of the spin depend on the history 
of the fluctuating field through the integrals of correlation functions. But the decay rate for the white noise case does 
not depend on history. Based on this observation, we say that spin system has memory if the fluctuating field driving 
it has finite correlation time. This definition is consistent with a rigorous notion of memory, since S(t), understood as 
a stochastic process, is Markovian (i-e. memoryless) if and only if the noise driving it is white noise. Memory in 
the system has important consequences, because it allows control of spin decay, through the application of external 
field. 

The role of memory can be understood as follows. Consider first the case of no constant field. Each realization 
of the fluctuating field uj(t) determines corresponding realization of (S(t)) = U(t,0)(S) via the equation of motion, 
where co(t) defines an instantaneous rotation axis for the spin vector. Hence, each realization draws a trajectory of 
the spin on a sphere of radius \{S)\. The evolution can be viewed in a simplified "stroboscopic" picture in which, 
within a correlation time r c , the axis co(t) does not change much, and subsequently, it jumps to a new random value 
within a range on the order of loq. Eventually such a sequence of rotations around randomly selected axes "smears" 
the trajectories of the spin over the whole sphere. This results in the decay of the spin averaged over all realizations. 

The situation is different when a constant field along, say, the z-axis is present. This becomes clear when we view 
the evolution in the reference frame rotating with frequency fio around the z-axis. Consider a particular realization 
of the fluctuating field uj(t), and denote it in the rotating frame as uj'(t). Since ui' z (t) = ui z (i), the fluctuations in the 
direction parallel to the constant field are unchanged when shifted to the rotating frame. On the other hand, the 
projection of the fluctuating field on the plane perpendicular to z-axis, = (uj x (t),u)y(t),0), is modified when 

viewed in the rotating frame. Within a time of order the correlation time r c , the projection uj'j_(t) revolves around the 
z-axis with frequency and if £Iqt c 3> 1 , 3'±[t) is able to perform many revolutions between jumps. The effect of 
the jumps becomes negligible in comparison with the systematic rotation. Thus the stochastic character of dj' ± (t) is 
lost, and the only contribution to the decay caused by fluctuations comes from oj z (t). A similar effect provides stable 
motion of a spinning top, and for that reason we will refer to it as a gyroscopic effect. If Qqt c <C 1, then revolution in 
between jumps is negligible, so ijj_(t) ~ u)±(t), and by necessity the decay induced by the fluctuations must be almost 
the same as if there was no constant field at all. 

The suppression mechanism of the decay described above is inoperative in the white noise limit where the correlation 
time vanishes, so the fluctuations of the rotation axis are instantaneous. Hence, the decay caused by white noise 
cannot be affected by a constant field of finite intensity. Moreover, even if the strength of the fluctuations is small in 
comparison to Q,q, the effect of fluctuations cannot be neglected if the correlation time is too short. 

Note that the important dimcnsionless parameter here is flo T c, and the conclusion of this paragraph remain valid 
even for non-Gaussian noise as long as it has nonvanishing correlation time. We verify this assessment in the following 
section. 

VI. NON-GAUSSIAN PROCESS: TELEGRAPH NOISE 

In this section we present an example of a non-Gaussian process, telegraph noise, i.e., dichotomic noise corresponding 
to a memoryless continuous-time stochastic process that jumps between two distinct values, sometimes called burst 
noise [l8|,[33|. We consider telegraph noise in one component of the field, ui z (t). The Rabi frequency of the fluctuating 
field, w z (t), jumps randomly between +loq and — ojq, with a jump rate w = (2t c ) . The average value and the 
correlation function of such a process is given by 

w^Z) = 0, Lo z {t)u z {t>) = c^e" 21 "!*-*'! =^e"^. (37) 

t-t' 

Because telegraph noise is non-Gaussian, the correlation function n(t — t') = to^e T " does not contain all the 
information regarding the process. 

Telegraph noise can model the effect of a spin 1/2 impurity atom near the spin [3, S @- The impurity spin has 
probabilities p + and p_ for being in the spin-up and spin-down state, and a hopping rate for transferring from one 
spin state to the other. Neglecting the back-action of the magnetic moment on the impurity, this model reduces to 
telegraph noise. 

Our goal is to find the average evolution matrix U (t, 0) for the telegraph noise case. The technique we have used 
previously to find the average does not work in the case of telegraph noise, because unlike with Gaussian processes, 
the closed form of the cumulant generating functional is not known. Fortunately, it is still possible to find an exact, 

analytical formula for U(t, 0) if the fluctuations are parallel to the constant field. One way of solving for U(t, 0) is 
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to follow Ref. [3J|. We can write the average evolution operator as U(t,0) = U(t,0), + U(t,0)_, where U(t,0) ± are 
the averages conditional on ui z {t) being equal to ±ujq. These averages vary in time because of equation of motion (J7]), 
and because uj z (t) jumps. Thus, 



-u(t,o) + 



(It 
d 



(n o +uJo)i z U(t,0)_ 



U(t,0)_ 



1 
1 

2^ 



U(t,0) + + —U(t,Q)_, 



U(t,0)_ 



1 

2^ 
1 

2^ 



J7(t,0). 



Adding and subtracting these equations, we obtain 



-C/(i,0) = Sl i*U(t,0)+u> o i z AU(t,0), 



—AU(t,0) = n Q e z AU{t,0) +uj o i z U(t,0) AU{t,0), 

dt T r 



(38) 
(39) 

(40) 
(41) 



where AU(t, 0) = U(t, 0) + — U(t, 0)_ . The next step is to eliminate AU(t, 0) from Eq. (|4*D| by substituting the solution 
to Eq. (|4"Tj) . with the initial condition J7(0,0) ± = 1/2. This yields a closed equation for average evolution matrix: 



dt j 



U(t,0) = n o e z U{t,0) + uj 2 I dt' e- t,/T "(i z ) 2 U(t-t\0) 



(42) 



Note the particular form of the "memory integral" in the equation (|42[) , where the dynamical variable is evaluated at 
times previous to t. In this case the system clearly has memory and it remembers not only the history of the noise 
but also that of its own. It is easy to check that the solution to this equation which satisfies the initial condition 

U(0, 0) = 1 is given by 



U(t,0) 



,-t/2r 



cosh ( 7— J +7 1 smh(j-^ r 



R(n t,z), 



(43) 



where 



V^4^ 



2 2 
T c 











7=1 V 1 - 4w o T c I ■ ( 44 ) 

V o 01/ 

The average expectation values of the spin are obtained by applying the average evolution matrix (|4"3")l to the initial 
vector of spin operators, (S(t)) = U(t,0)(S(0)), which yields, 



(S x (t)) 



(Sy(t)) 



cosh 



cosh 



1 - ^Irl 



t 
2r r 



sinh 



1 - Au 2 T 2 — 



s/i - ^ 



sinh ( y/1 - 4cjgr| 



2r, 



[cos{n Q t)(S x ) - sm(n t)(S y )} , (45) 



[cos{n t)(S y ) + sm(n t)(S x )} , (46) 



(S z (t)) = (s z ). 



(47) 



Neither constant field Oo; nor field fluctuations ui z (t) couple to the z component of the spin, hence it remains constant 
throughout the evolution. The effects of the presence of £7 an d w z (t) on the components of the spin in xy plane 
are simply superposed: constant field causes trivial precession around z axis, while the noise makes the x and y spin 
components decay on average. Moreover, if the constant field and the fluctuating field are parallel, the gyroscopic 
effect can not be observed. The only memory effect present in this case is a quadratic behavior in t of the decay rate 
at t <C r c . Figure [5] plots (S x {t)) versus time for telegraph noise in ui z (t), for the case of no static external magnetic 
present. As w t c increases, (S x (t)) decays more quickly, but for ujqt c > 1/2, the square root yl — 4ojqT^ becomes 
imaginary and the solution changes form from pure decay to decay with oscillations. 

It is of interest to see what happens when noise is present in the component of the field perpendicular to the constant 
field f2oj e.g., in uj x (t). Unfortunately, in this case, it is not possible to obtain an analytical formula for the average 
evolution matrix. According to the previous discussion, one would expect that, since telegraph noise is colored noise 
(the correlation time r c is finite), the gyroscopic effect should work in a fashion similar to the Gaussian colored noise 
example. This is indeed the case. Numerical solutions of (S z (t)) versus time are plotted in Fig. [5] and suppression of 
the decay is more efficient with increasing Qqt c . 
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-0.4; 



2 4 6 8 10 12 14 

tlT c 



FIG. 5: {S x (t)) versus time for a spin subject to telegraph noise in only the z-component of the field, uJ z (t), and £1q = 0. (S x (t)) 
is identical. Nonvanishing values of flo simply cause precession of the spin about the z axis. 




FIG. 6: (S z (t)) versus time for a spin subject to telegraph noise only in the x direction, with ojqt c = 1, perpendicular to the 
constant field S7o which is in the z direction. As the parameter f2or c is increased (see legend), suppression of the decay of the 
spin increases. 



VII. SUMMARY AND CONCLUSION 



We introduced a stochastic model for spin in the presence of a deterministic magnetic field and a fluctuating magnetic 
field (noise), and we derived the stochastic equations of motion for the spin vector (S(t)). The environmental degrees 
of freedom, represented by the fluctuations of the field, are eliminated by taking the average over the fluctuations. 
We made the external noise assumption, wherein no back-action of the system on the environment is present. The 
most pronounced consequence of this assumption is that the system does not go into equilibrium with a thermal 
environment, but instead goes to the most democratic density matrix, hence the system tends towards steady state of 

the completely mixed state, and (S(t)) — >• as t — > 00. This is in contradistinction to the case of noise experienced by 
a system in contact with a bath, wherein the bath affects the system and the vice versa. The lack of mutual interaction 
(back-action) means that the fluctuation-dissipation theorem [ItJ cannot be applied, and a thermal equilibrium state 
of the system cannot be achieved at large times. If back-action on the bath due to the system is present, the source 
of the fluctuations and the dissipation is the same, they are connected through fluctuation-dissipation theorem, and 
thermal equilibrium of the system must result. 

We explicitly considered Gaussian colored noise, Gaussian white noise and non-Guassian telegraph noise. From 
our studies we conclude that if the system has memory, i.e., it has a finite correlation time, then the system can 
by manipulated by means of an external magnetic field, and the dccohcrcncc induced by the fluctuations can be 
significantly suppressed. In contradistinction, in the white noise limit, when the correlation function tends to a Dirac 
delta function, the decoherence cannot be controlled in this way. We introduced simple and intuitive considerations 
for why the suppression of decay can be achieved for any stochastic process as long as the correlation time is long 
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enough. 

In futurepublications we will present a quantum master equation approach based upon the techniques developed in 
this work [23| , wherein quantum master equations are derived an yield an equation for the density matrix for the spin. 
Before averaging over the random magnetic field, the equations of motion are stochastic because of the coupling of the 
stochastic magnetic field to the spin. Averaging over the stochastic processes can be calculated analytically to yield a 
deterministic equation of motion for the density matrix averaged over the stochastic magnetic field realizations. We 
will also discuss the effects of stochastic fields on two-qubit systems that are and are not entangled, to see how the 
correlation and entanglement of the qubits are affected by different kinds of noise [36] . Moreover, we will explicitly 
explore the effects of magnetic noise on magnetometers and atomic clocks 37]. 
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